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Abstract. — We analyze the convergence of the exponential Lie and exponential Strang splitting 
applied to inhomogeneous second-order parabolic equations with Dirichlet boundary conditions. A 
recent result on the smoothing properties of these methods allows us to prove sharp convergence 
results in the case of homogeneous Dirichlet boundary conditions. When no source term is present 
and natural regularity assumptions are imposed on the initial value, we show full-order convergence 
of both methods. For inhomogeneous equations, we prove full-order convergence for the exponential 
Lie splitting, whereas order reduction to 1.25 — e for arbitrary small e > for the exponential Strang 
splitting. Furthermore, we give sufficient conditions on the inhomogeneity for full-order convergence 
of both methods. Moreover our theoretical convergence results explain the severe order reduction 
to 0.25 — e for every e > of splitting methods applied to problems involving inhomogeneous 
Dirichlet boundary conditions. We include numerical experiments to underline the sharpness of our 
theoretical convergence results. 



1. Introduction 

The focus of this work lies in a rigorous error analysis of the exponential Lie and Strang 
splitting applied to second-order inhomogeneous parabolic evolution equations 

d t w(t,x,y) = £{d x ,d y )w(t,x,y) + ij{t,x,y), (x, y) G ft = (0, l) 2 , t e]0,T], 

(1) w(0,x,y) = w (x,y), 

w(t,-,-)\ dQ = for all t € [0,T], 

where £(d x ,d y ) = d x (a(x,y)d x ) + d y (b(x,y)d y ) denotes a strongly elliptic differential operator 
with sufficiently smooth (positive) coefficients a and b. In particular, we elucidate the following 
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phenomena: for arbitrary inhomogeneities tp the exponential Strang splitting naturally induced 
by the decomposition of C suffers from an order reduction whereas the exponential Lie splitting 
is full-order convergent, see Figure [T] in Section 5 for a numerical example. Nevertheless, for a 
certain choice of ip also the Strang splitting reaches full second-order convergence, see Figure |3| 
A recent result on the smoothing properties of exponential splitting methods, see [12J , allows us 
to prove this order behavior rigorously. We can show full-order convergence of the exponential 
Lie splitting method, despite the fact that its local order of consistency reduces to 0.25 — e for 
arbitrary small e > 0. Furthermore we can state a sharp result on the global order reduction to 
1.25 — e for arbitrary small e > of the exponential Strang splitting applied to equations with 
general inhomogeneities ip. 

Operator splitting methods have gained a lot of attention in recent years, as in many situations 
they allow a keen reduction of the computational cost. For a general introduction to splitting 
methods, we refer to [7J, [llj and references therein. Particularly, exponential splitting methods 
with relatively bounded operators were studied in [8 J ; the situation of two unbounded operators 
in the homogeneous parabolic case was studied in [4] and [5]. However, due to the recent 
smoothing result, we can improve the error bounds in the homogeneous parabolic setting given 
in the last two references. Resolvent splitting applied to inhomogeneous equations was discussed 
in [T3] , 

This work in particular distinguishes from the above as inhomogeneous Dirichlet boundary 
conditions are included in the analysis. It is well know that inhomogeneous Dirichlet bound- 
ary conditions generally lead to a severe order reduction of splitting methods^ where for an 
introductory reading we refer to [7J Section IV. 2]. Our theoretical analysis allows us a pre- 
cise explanation of these phenomena, and numerical experiments confirm the sharpness of the 
derived convergence bounds. 

For the analysis, the model problem ([!]) is formulated as an abstract evolution equation in 
L 2 (Q), where for convenience we denote by || • || the L 2 (0) as well as the L 2 (0)-operator norm. 
Thus, the model problem reads 

u'{t) = Lu(t) + g(t) = Au(t) + Bu(t) + g(t), t e]0, T], 

u(0) = u 

with u(t) = w(t, •, •), Au = d x (ad x uj, Bu = d y (bd y u\ and g(t) = ip(t, -, •). 

For the numerical time integration, we apply the so called dimension splitting, where the 
differential operator L is split along its dimensions. Let h > denote the step size. At time 
tn+i = (n + l)h the exact solution of ^ is approximated by exponential Lie splitting 

(3) u n+l =e hA e hB {u n + hg(t n )) 
or the exponential Strang splitting 

(4) u n+1 = M£ B (£ B Mu n + hg(t n + I)) , 



1 Particular attention will be paid to the order reduction caused by boundary conditions since that is often the main 
reason for a disappointing convergence behavior with splitting methods; [7J p. 348]. 
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respectively. In Section [3] the advantage of choosing this particular structure of the exponential 
Strang splitting is discussed in more detail. 

In the bounded setting one can easily see via Taylor series expansion that the classical order 
of the exponential Lie splitting ^ and the exponential Strang splitting Q is one and two, 
respectively. Under particular assumptions on the initial value uq and inhomogeneity g we will 
show that full-order convergence of both methods also holds true in our unbounded setting, see 
Section [3j Theorems [6] and [7J respectively. In particular we can explain the numerical observed 
full-order convergence of the exponential Strang splitting in [5] . In the situation of an arbitrary 
inhomogeneity g, we will prove full-order convergence for the exponential Lie splitting, whereas 
order reduction to 1.25 — e for arbitrary small e > for the exponential Strang splitting, see 
Section[3j Theorems [6] and [7j respectively. In Section[4]we explain the disappointing performance 
of splitting methods applied to equations involving inhomogeneous Dirichlet boundary conditions 



by applying our newly derived theoretical results, see Corollaries 10 and 11 Finally, in Section 



[5] we illustrate the sharpness of our theoretical convergence results with numerical experiments. 
We commence with a section on the analytic framework. 



2. Analytic framework 

In order to start with the convergence analysis we need some functional analytic results on 
the domains of the operators in ([2]). In the following we set 

v(L) = H 2 (n)nH^{n), 

V{A) = {ue L 2 {n); d xx u,d x ue L 2 (n), u(0,y) = u(l,y) = f.a.e y G (0,1)}, 

V{B) = {u G L 2 (n) ; dyyu, d y u G L 2 (n), u(x, 0) = u(x, 1) = f.a.e x G (0, 1)}. 

We will also need the compositions of the operator L and therefore introduce the spaces 

V{L k ) = {u G V(L k ~ x ) ; Lu G D(L fc_1 )}, k = 1, 2, 3, . . . 

Note that all the appearing operators (V(A),A), (V(B),B) and (V(L),L) generate analytic 
semigroups of contractions on L 2 (il), see |12j . For a general introduction to the theory of 
analytic semigroups it is referred to 1 14! . Furthermore note that D(L 7 ) is free from boundary 
conditions if and only if 7 < \ , see [2] . 

In the following we state some results on the regularity and compatibility between the domains 
of the full operator L, the split operators A, B, and their compositions, which we will need later 
on in our convergence analysis. 

Lemma 1. — The following regularity results hold 

V(L 2 ) C H A (Q) and V(L 3 ) C H 5 - £ (Q) for all e > 0. 

Proof. — This result can be derived from an analysis of the classical compatibility conditions 
for a strongly elliptic operator with Dirichlet boundary conditions on a square domain, see for 
instance [3]. A complete proof can be found in [6]. □ 
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As the operator (T>(L),L) is considered on the unit square f2, i.e., a corner domain, the 
standard regularity results do not hold true. This is due to the so called corner singularities 
that arise when solving partial differential equations on non-smooth domains. In particular, for 
general coefficients, T>(L 3 ) £ H 6 (n). 

Lemma 2. — The domains of the operators satisfy the following compatibility condition 

V(L 2 ) C V(A 2 ) n V(AB) n V{BA) n V(B 2 ). 

Proof. — From Lemma [T] we have 

V(L 2 ) C {u G H\n) n Hq(U) ■ Lu\ 9n = 0} . 

Let u G V(L 2 ). Then Au, Bu G H 2 (Q). As Q C M 2 , we have H 4 (n) n flj(n) C C 2 (H) n ^(fi) 
using standard Sobolev embeddings results. Hence as u G H 4 (Q) n i?o(^) we nave f° r ^o>yo G 
{0,1} 

^4w(sc, 2/)| y=yo = d x (a(x,y )d x )u(x,yo) =0, 

Sm(x, y) U=x = <^ (6(z , y)9 y )u(x , y) = 0. 

Thus, as = Lit|gn = Au\qq + i3u|ao we have by using ([5]) ^4u|an = Bu\qq = 0, i.e., Au,Bu G 
P(i)nD(S). " " □ 

Lemma 3. — Let f G £>(L 2 ) C I>(AB) n V(BA). Then for every A > the solution u to the 
elliptic problem 

(XI - A)u = f in O, 
^ ^ u(0, y) = u(l,y) = for almost every y G (0, 1) 

lies in V(AB) DV(BA). 

Proof — For u solving @ and / G £>(L 2 ) C we have u G see [12]. Thus, 

^lu = Au — / G T>(B), i.e., u G T>(BA). Next we show that u G T>(AB). For this, we compute 

(7) 5« = fy(ftfyu) = dybd y ((XI - A)" 1 /) + ^((AI - A)- 1 /). 

The relation d y ( XI- A) ^ = (XI - A)~ l d x ((d y a)d x ) (XI - A)~ l (see [H Section 7]) shows that 

d y (XI - A)-\ d yy (XI - Ay 1 ■ L 2 (Q) -> V(A) c L 2 (0) 

for smooth coefficients a. Consequently Bu G T>(A), which proves the Lemma. □ 

The above Lemma in particular implies that for all A > 

(XI - A)- 1 ^^, (XI - B)- 1 ^^ : V(L 2 ) C L 2 (tt) -> V(AB) n V(BA) C L 2 (Q). 

Thus, by using Cauchy's integral formula we can conclude that 

etA \v(L*y etB \v(L*) : v ( l2 ) c ~> n c 

uniformly in i > 0. Hence, by the closed graph theorem there exists a constant C such that 
for any h > 0, sup a>T6[0jh ] ||L4, B]e aA e rB L _2 || < C. Furthermore, as for < 7 < |, V(L^) 
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is an invariant subspace under the semigroups generated by the split operators, we have for 
< e < \ that sup SiTe r 0j/l i \\L~ l+£ [A, B]e sA e TB L~ l ~ e \\ < C uniformly in h > (for some 
constant depending on e). 

Next we study the regularity of the solution to our model problem Q. Its exact solution is 
given by the variation-of-constants formula 

u(t) = e tL u + [ e {t - r)L g(T) dr, < t < T. 
Jo 

At time t n+ \ = t n + h, with a step size h > 0, the solution can be rewritten as 



rh 

i(t n+1 ) = e hL u{t n ) + / e {h - s)L g(t n + s) ds. 
Jo 



Taylor expansion of g(t n + s) at t n yields 

rh . rtn+s . 

u(t n+1 ) = e hL u(t n ) + e {h - s)L (g(t n ) + sg'(t n ) + jf (t n + s - t) 9 "{t) drjds. 

We introduce the following notation (see [3]): 

Definition 4- — For h > we define the bounded operators Xj by setting 

A-- — f h e {h - s)L ds 7>1 

J ~yJ e (j-l)! dS ' 3 ~ l 

and Ao = e hL . We likewise define Qo = e hA , /3q = e hB and for j > 1 

cu = i [ h e^ A 7 S —^ds, fa = i f k ^ B 7 S —-ds. 

3 h> J u-iy. ^ wy 

With this notation the exact solution to ([2]) can be written as 

(8) m(Wi) = A n(t n ) + h\ ig (t n ) + h 2 X 2 g'{t n ) + E 3 ,n(M> 
with the remainder 

fft / f-tn+S 

(9) i? 3 ,n(/l) 



Jo e(h ~ S)L Ut {tn + S ~ t)9 " {t) dT ) ds ' 



We note for later use that following relations hold on L (O) 

(10) Xj = -I + hLX j+ i, j>0, 

which can easily be seen by integration by parts. 

We now state some well known regularity results for analytic semigroups, see |10j for further 
details. This allows us to play back the regularity assumptions on the exact solution to regularity 
assumptions on the initial data. For some < 7 < 1 and 9 > let the data of problem ([2]) 
satisfy 

(11) uoeV(L^), Lu + g(0)eV(L^), g G C 1+e ([0, T}; L 2 (n)), 
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then the solution of ([2]) possesses the regularity 

(12) ueC 2 ([0,T];L 2 (n))nC\[0,T};V(L^)). 

Note that if we consider second-order methods, we need more regularity of the solution. For 
analytic semigroups, this requirement can again be easily expressed in terms of the data. If, for 
some < 7 < 1 and some 9 > 0, the data satisfies 

9 eC 2+fl ([0,T];L 2 (O)), 

u Q £V(L~>), Lu + g(0) GD(LT), L 2 uq + Lg(0) + </(0) G 2?(L 7 ), 
then the exact solution of the evolution equation Q satisfies 
(14) u G C 3 ([0, T]; L 2 (Q)) n C 2 ([0, T]; D(L 7 )). 



3. Convergence results 

In order to state our sharp convergence results we need the following lemma about the smooth- 
ing property of splitting methods recently derived in [12J. 



Lemma 5. — For < a < 1, the following smoothing properties hold 

\\(-L) a (e hA e hB ) n \\ < Ct~ a , < t n = nh < T, 

\\(-L) a (e% A e hB e% A ) n \\ < CC\ < tn = nh < T, 

where the constant C can be chosen uniformly on [0, T] and, in particular, independently of 
n > 1 and h > 0. 

Now we are able to state the main convergence results for the exponential Lie and Strang 
splitting, see Sections 3J_ and 3^, respectively. 

3.1. The exponential Lie splitting. — As a numerical scheme, we now consider the expo- 
nential Lie splitting method for inhomogeneous equations, i.e., the numerical approximation to 
the exact solution at time t n+ \ = t n + h is given by 

(15) u n+ i = e hA e hB («„ + hg{t n )) , 

and initial condition uq = u(0). Under natural regularity assumptions in time and space on the 
inhomogeneity g we obtain first-order convergence. 



Theorem 6. — The global error of the exponential Lie splitting method (15) satisfies the fol- 
lowing bound: For all fixed T and 7 > ; there exists a constant C such that for all h > 0, we 
have 

(16) \\u(t n ) - un\\ < Ch sup (||LV(t)|| + \\L^g(t)\\ + \\g' 

te[o,T] v 

for all n such that t n < T, where u(t) is the exact solution of 
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We recall that the domain P(L 7 ) is free of boundary conditions for 7 < Thus full-order 
convergence of the exponential Lie splitting only requires some additional smoothness in space 
of the inhomogeneity g. 

Also note that if for some < 7 < 1 and 6 > 0, u G V(L 1+r ), g(0) € £>(L 7 ) and g £ 
C 1+e ([0, T]; L 2 {yi)), then the exact solution of problem ^ will satisfy the regularity assumption 
( p~2| ) , in particular sup tg [ T ] ||L 7 u'(i)|| < C. 



Proof. — In the following we will denote by S the exponential Lie splitting, i.e., S 
The exact solution of ([2]) is given by 

(17) u(t n+1 ) = X u(t n ) + h\ig(t n ) + R2, n (h), 

where we have set 

ft n +S 



^hA^hB 



(18) 



R2,n(h) 



Jo 



,(h-s)L 



g'(r) dr I ds 



In order to prove convergence of the method, we first take a look at the error e n +i = u rt +i 
u(t n+ i). Subtracting (17) from (15) leads to 

R2,n(h)- 



e n+ i = Se n + (S - X )u{t n ) + h(S - Xi)g{t ri 

I), as well as the equation in the form 



(19) 



Using the recurrence relation (10), namely Ai 



u{t n ) + L- l g{t n ) = L^u'it n ), we obtain 

(20) e n+1 = Se n + S (u(t n ) + hg(t n )) - \ L~ 1 g(t n ) + L~ 1 g(t n ) - R 2 , n ( h )- 

The convergence analysis of the exponential Lie splitting for the homogeneous problem stated 
in [4J proof of Thm. 3.4] allows us to express the identity operator / restricted to the domain 
V(L) as 



(21) 



\V(L) 



S(I -hL) + h 2 K 2 , 



where the operator K,2 is such that h^^L 1 || < C for some constant independent of h. Inserting 
the representation (21) of the identity operator in front of the terms with L~ l in (20) yields 

e n+ i =Se n + S (u(t n ) + hg(t n )) - I\ L~ 1 g(t n ) + IL~ 1 g(t n ) - R2, n (h) 

=Se n + S(I - A + hL\v)L- l u'{t n ) - h^XoL^u' \t n ) + /i 2 /C 2 L- 1 5 (t„) - R 2 , n (h). 



Using once more the recurrence relation (10) we obtain 

(22) e n+1 = Se n + /i 2 5(Ai - X 2 )Lu'(t n ) - h^X^u' \t n ) + /i 2 /C 2 L- 1 < 7 (t n ) - R 2 , n {h). 
Solving the above error recursion we obtain with the aid of Lemma [5] for < 7 < 1 

n— 1 n—k 

K+ill < h 2 WS^L^WWL-^JkW + h 2 \\J n \\ + \\S n - k \\\\R2,k( h )\\ 
h 



(23) 



fe=0 
n-l 



k=0 
n—k 



< 



k=0 



((n — k)hy-~< 



\L- l ^J k \\ + h 2 \\J n \\ + Y,\\R2,k{h)l 



k=0 
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where we have set 

J k = 5(Ai - X 2 )Lu'(t k ) - K 2 \ Q L- 1 u'{t k ) + K 2 L- X g(t k ). 
Using the domain invariance of 2?(L 1_7 ) for < 7 < 1 under the semigroups, we have 

\\L- 1+ -<J k \\ < c(||LV(t fc )|| + \\L-y g (t k )\\y 

Note that h\\K2L~ || is bounded and thus /j 2 ||J n || < Ch. Furthermore 
C/t 2 sup t6[0jT] \\g' 



R 



2,k 



< 



see (18). Thus, by estimating the sum in (23) by its Riemann inte- 

□ 



gral, we obtain the desired convergence result (16). 

3.2. The exponential Strang splitting. — In this section we are interested in applying 
a second-order exponential splitting method to the inhomogeneous problem ([2]). As numerical 
approximation u n +i to the exact solution u(t n +i) we choose the following scheme 



(24) 



u 



n+l 



e 2 A e 2 B l e 2 B e 2 A u n + hg(t n + 



Clearly there are many possible choices of applying second-order exponential splitting methods 
to inhomogeneous problems, such as for instant the scheme 



(25) 



U n +1 



e^ A e hB e-2 A 



u n + ^g(t n j) + ~g(t n + h), 



which we will call Strang B splitting in the following. The main motivation behind our choice 
( |24[ ) lies in the crucial fact that for uq G T^{L) the numerical solution u ra +i remains in the domain 
T>(L) for every n > 0, which can be seen from the estimate 

\\Lu n+ i\\ < \\L^ A e hB e'i A L- l \\\\Lu n \\ + \\hLe'i A ^ B \\ \\g[t n + |) || 

< C^||Lu n || + \\g(t n + ^)| 

that uses the domain invariance of T>(L) under the semigroups. In contrast, the numerical 
solution of Strang B (25) remains in the domain T>{L) only under the unnatural assumption 
that g(t) 6 T>(L) for all t > 0. The favorable choice of (24) is also confirmed in the numerical 
experiments given in Section [5] 



Theorem 7. — The global error of the exponential Strang splitting method (24) satisfies the 
following bound: For all fixed T, 7 > 0, and v with < v < 1, there exists a constant C such 
that for all h > 0, we have 



(26) 



\u{t n ) 



'11 r. 



< Ch 



1+1/ 



sup 

te[o,T\ 



\L~<u"(t)\\ + \\L^g(t)\\ + \\^g'(t)\\ + \\g" 



for all n such that t n < T, where u(t) is the exact solution of 

Remark 8. — Note that full-order convergence of the exponential Strang splitting (24) requires 
the additional assumption that g G C([0, T]; P(L 1+7 )), for some 7 > 0. In particular the 
inhomogeneity ip(t, •, •) = g(t) needs to vanish at the boundary, i.e., tp(t, •, -)\dfi = 0. 
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Also note that if for some < 7 < 1 and 6 > 0, u £ V(L 2+ ~<), g(0) £ V(L 1+ ~<), g'(0) G V(L^) 
and g G C 2+£, ([0, T]; L 2 (J7)), then the exact solution of problem ([2]) will satisfy the regularity 
assumption ([14]), in particular sup tg [ T ] ||W(t)|| < c. 



Remark 9. — Full-order convergence of the exponential Strang splitting (24) without unnat- 
ural assumptions on the inhomogeneity g(t) can be obtained when considering the evolution 
equation ([2]) in the extrapolation space (T>(L)*, WL^ 1 • ||). We refer to pQ for an introduction to 
scales of Banach spaces. Precisely, we can adapt the previous proof and show that the global 
error measured in the \\L~ l • || norm satisfies the following bound for any 7 > 0: 

\\L-\u(t n ) - OH < Ch 2 sup (||LV(t)|| + \\L-y g (t)\\ + \\L-r- l g'(t)\\ + [|L-V'(*)||) 

te[o,T] v 7 

with a constant C that can be chosen uniformly with respect to h on [0,T]. 

Proof. — Let S = e^ A e hB e^ A . We start off with the analysis of the local error e„,+i = u n+ \ — 
u(t n +i). Subtracting the exact solution Q from the numerical approximation (24) we obtain 

e n+ i = Se n +(S- X )u{t n ) + h(e^ A e^ B - Xi)g(t n ) 

(27) 

+ h 2 {\^ A ^ B -\ 2 )g'{t n )-R,, n {h) 
with i?3 n given in @. In the following we set 

E n {h) = (S- X )u(t n ) + h(ev A ^ B - \ x )g{t n ) + h 2 ( l -e^ B - \ 2 )g'(t n ). 
With the identities Ai = ^(Ao — I)L~ l and A2 = (Ao — I — hV)L~ 2 we obtain 

E n (h) = Su(t n ) - X L- 2 u"(t n ) + h^ A ^ B g(t n ) + L^gitn) 
(28) u2 h , 

+ \£ A £ B g'{t n ) + L- 2 g'(t n ) + hL-'g'itn). 

Note that the following representations of the identity operator are valid on T>(L) and T>(L 2 ), 
respectively, see [4], 

I\v(L) = S(I-hL) + ^JC 2 , 

(29) 

I\ V(L2) = S(I - hL + ^L 2 ) + h 3 JC 3 , 
where the operators IC2 and /C3 satisfy the following bounds: for all v with < v < 1 

(30) h 2 \\K2L- 1 f\\<Ch 1+v \\L u f\\, 
and for all 7 > and v with < v < 1 

(31) h^WL^'K^fW < Ch 2+U \\L^ +U f\\ 
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for some constants independent of h and for all functions /. Using the first representation of 
the identity operator given in (29) in front of the terms with L _1 and L~ 2 in (28) we obtain 



E n (h) ={S - X )L- 2 u"{t n ) + h^ A e^ u {I - e^e^ A )g(t 



A IB i 



(32) 



+ ^K 2 L- l {g{t n ) + L-V(i n ) + hg'(t n )) + f el A ef B g'(t n ) - h 2 Sg'(t n ). 



Prom the convergence analysis of the exponential Lie splitting with h replaced by - we know 
that 

(33) I - Q2 a e 2 A = --ez^e? L + — fC 2 , 



see (21). Furthermore, by [3] we have for all functions / and all v with < v < 1 
(34) \\(S-\ )L- 2 f\\<Ch 2+v \\L v f\\ 



for some constant C independent of h and /. Thus, using (33) as well as the bounds (30) and 



( [34) ) in (|32j), we obtain that for < v < 1 

(35) \\E n (h)\\ < Ch 1+V (\\u"(t n )\\ + ||^(*n)|| + ||</(Ull)- 

Using the second representation of the identity operator given in ( 29 ) in front of the terms with 
L^ 1 and L~ 2 in (|28l) we obtain 



(36) 



E n (h) = (S- \ )L- 2 u"(t n ) + he h * A £ B (l + £ B £ A {\h - I))g{t n ) 

+ h ^ A ^ B (/ - et^et A )g'(t n ) + h^L^g{t n ) + h^L- 2 g\t n ) 
+ h 4 JC 3 L- 1 g'(t n ). 



From the convergence analysis of the exponential Lie splitting with h replaced by - we know 
that 



(37) 



I + ^B e \A(h L _ I 



/c 2 



I — e2 e2 



^ B e^ A L+^JC 2 , 



see again (21). Thus, inserting the relation (|37|) into (|36|) we obtain 

-2// 



(38) 



L\,(h) - (S - X )L- 2 u"(t n ) + ^ A ^ B JC 2 g(t n ) - ^SLg'(t n ) 
+ ^e h 2 A e~2 B IC 2 g'(t n ) + h^L- 1 g{t n ) + h 3 JC 3 L~ 2 g' (t n ) 



+ h^L- l g'{t n ). 

From the domain invariance of T>(L) under S and [5], we know that for < 7 < 1 

\\L- 1+ ~<(S - X )L~ 2 u"(t n )\\ < Ch 3 \\L^v!'{t n )\\. 
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Thus, using the domain invariance of the domain T>(L l 7 ), < 7 < 1 under the semigroups 
(38), as well as the bound given in (31), we obtain for < v < 1 the following estimate 



m 



(39) 



\L- 1+ ~<E n {h)\\ < Ch 2+ »(\\Vu"{t n )\\ + W +V g(t n )\\ + Wg'{t n ) 



Thus, using the bounds given in ( [35]) and (39) as well as the smoothing property for the expo- 
nential Strang splitting, see Lemma |5| we obtain for < v < 1 



l e n+l| 



n— 1 n— k 

< £ H^-^-iiiL-^fc^ii + \\E n (h)\\ + £ \\s n 

k=0 k=0 



\\R3,k(h)\\ 



n-1 



(40) 



^ Chl+U E ((n _fc )/l) i- 7 (ll LV/ fa)H + W L ~ <+U 9(tk)\\ + ||LVftb) 



fc=0 



n—k 

+ Ch 1+V (\\u"{t n )\\ + Wg(t n )\\ + || 5 '(i n )||) + Y, 

k=0 

By ([9]) we have ||i?3 ) fc(/i)|| < C/i 3 sup tg [ 0j T] Hence, approximating the sum in (40) by its 
Riemann integral, leads to the desired estimate ( |26[ ). □ 



4. Application to inhomogeneous Dirichlet boundary conditions 

The main aim of this section is to understand the disappointing behavior of splitting methods 
when applied to problems involving inhomogeneous boundary conditions with the aid of the 
convergence results derived in the previous section. Our model problem in this section thus 
reads 

d t w(t,x,y) = C(d x ,d y )w(t,x,y), (x, y) e fl = (0, l) 2 , t e]0,T], 
(41) w(0,x,y) = w (x,y), 

w (V,0| 9n = /(*,v) for all te [0,T], 

where C(d x , d y ) = d x (a(x, y)d x )+d y (b(x, y)d y ). In the following we denote by Gf(t) the Dirichlet 
map, i.e., the solution to the elliptic problem 

C(d x , d y )v = in Q, 
v\an = f(t). 



In Section 4.1 we will discuss the case where the Dirichlet map Gf(t) is analytically known and 
smooth in space and time, whereas in Section 4.2 we focus on the more common situation where 
/ is only known as a boundary function. 
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4.1. Formulation as an inhomogeneous equation in L 2 (0). — Here we assume that the 
Dirichlet map Gf(t) =: F(t) is analytically known and in particular smooth in space and time, 
i.e., our model problem (41) satisfies 



(42) 



W(t,; 



F(t,-,-)\ m for all t G [0, T\. 



The ansatz is to rewrite (41), (42) as an inhomogeneous evolution equation in L (ft) 
with homogeneous Dirichlet boundary conditions. Employing the transformation U(t, x, y) = 
w(t, x, y) — F(t, x, y), we see that U is the solution to the problem 

dtU(t,x,y) = C(d x ,d y )U(t,x,y) + i(>(t,x,y), (x, y) G fi = (0, l) 2 , t e]0,T], 
(43) U(0,x,y) = U (x,y), 
U(t,;-)\ m = 

with ip(t,x,y) = C(d x ,d y )F(t,x,y) - d t F(t,x,y) and U (x,y) = w (x,y) - F(0,x,y). This 
equation can be formulated as an abstract evolution equation in L 2 {Q) with u(t) = U (t, •, ■) and 
g (t) = ^(t,-,-), i.e., 

u'(t) = Lu(t)+g(t), te]0,T], 
u(0) = uq, 



(44) 



where now the problem is equipped with homogeneous Dirichlet boundary conditions, i.e., 
T)(L) = H 2 (£l) n Hq(Q). Thus, we can apply the convergence results for the exponential Lie 
and Strang splitting stated in Section [3} 

We are in particular interested in how the regularity assumptions on the initial value uq and 
the inhomogeneity g of the reformulated inhomogeneous problem (44) relate to the regularity 
assumptions on wq and the boundary data F of the original inhomogeneous boundary value 
problem (41). Therefore we introduce the operator L F ^ = C{d x ,d y ) with the domain 

V{L m ) = G H 2 (n) ; (p\ an = F(0, •, -)\ dn } . 

First-order convergence of the exponential Lie splitting requires the following regularity as- 
sumption on the initial value uq and inhomogeneity g 

(45) noGP(L^), 1/(0) = Lu + g(0) eV(L^), g G C 1+e ([0, T\; L 2 (n)) 

for some 7 > 0, see Theorem [6j Expressed in terms of the original data in (41), assumption (45) 

I* GclCIS 

u =w -F(0) GV(L^), 

u '(0) = w '(0) - F'(0) = L F{0) w(0) - F'(0) G V(L~<), 
C(d x , d y )F(t, ■, •) - d t F(t, ■, •) G C 1+e ([0, T];L 2 (£l)). 

Thus, for first order convergence of the exponential Lie splitting method, the initial value of the 
original problem (41) needs to satisfy for some 7 > 

(46) w G V(L m ) n H 2+2 ^(Q) = {<pe H 2+2 \U) ; plan = F(0, •, -)|an} . 



EXPONENTIAL SPLITTING METHODS 



13 



Hence upon using Theorem [7j a convergence rate of order 1.25 — e of the exponential Strang 
splitting method can be achieved under the additional regularity assumptions 

gGC 2+e ([0,T];L\n)), 

n G V{p), u'(0) = Lu + g(0) G V{p), 

u"(0) = L 2 u + Lg(0) + ff'(O) G V(L"<) t 



for some 7 > 0. In terms of the solution w and boundary data F of the original problem (41) 
the additional assumptions read 

u"(0) = w"(0) - F"(0) = L 2 F{0) w(0) - F"(0) G V{p\ 

C(0 X , d y )F(t, ■, •) - d t F(t, ■, •) G C 2+e ([0, T]; L 2 (ft)). 

Thus, convergence of order 1.25 — e requires the additional smoothness assumption on the initial 
value wo 

(47) u> G P(4(o)) n # 4+27 (^) = G i7 4+2 ^(^) ; C(d x , d y )<p\ aQ = <p\ dn = F(0, •, -)\ aQ } 
for some 7 > 0. 

Note that in order to have full-order convergence of the Strang splitting applied to inhomo- 
geneous equations, the inhomogeneity g needs to lie in the domain 2?(L 1+7 ) for some 7 > 0, 



precisely g G C([0,T];"D(L 1+7 )), see Remark [8j In terms of the original problem (41) this leads 
to the following assumption on the boundary data F 

C(d x ,d y )F(t, ■, •) - d t F(t, ; •) G C([0,T]; H 2+2 ~<{Sl) n H^O)) 

for some 7 > 0. We summarize our findings in the following corollary. 



Corollary 10. — The exponential Lie and Strang splitting applied to (42) are convergent of 
order 1, respectively, 1.25 — e (for every e > 0) under the natural assumptions (46) and (47), 
respectively. Full second-order convergence of the exponential Strang splitting method is obtained 
if the boundary data F satisfy 

d t F(t,;-)\ dn =£(d x ,d y )F(t,;-)\ m 

for all t G [0,T]. 

4.2. Formulation as an inhomogeneous equation in the dual space (D(L)*, • ||). 
— Here the Ansatz is to rewrite the original problem (41) as an abstract evolution equation 
in the extrapolation space TJ(L)* with the corresponding \\L-~ 1 ■ \\ norm. First we formulate the 
problem in L 2 (ft) as follows 

(48) L- 1 y'(t) = y(t) + Gf(t), 

y(0)=y = w(0), 

see [9] for further details. In order to analyze the convergence of the splitting methods applied 
to ( 48 ) , we need to reformulate the equation as an evolution equation in the extrapolation space 
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(P(L)MIL- 1 • ||), i.e., 

( . Q , y'(t) = Ly(t) + LGf(t), 

2/(0) = yo- 

Note that this is the usual approach if the Dirichlet map of the boundary data is not known 
analytically. For instance, in a simple one dimensional situation with constant coefficients, 
namely 

dtw(t,x) = d xx w(t,x), a; €(0,1), t e]0,T], 

w(Q,x) = wq(x), 

w(t,0) = a(t), w(t, 1) = P(i), 

employing standard second-order central differences with grid size Ax, yields a semidiscrete 
problem of type 

W'(t) = ^ diag [1, -2, l]W(t) + [a(t),0, . . . , 0, /3(t)] T , 



which is of the form ( 49 ) 
F 

(50) 



First, we analyze the convergence in the dual space (T>(L)* , \\L 1 • ||) and consider 

u'{t) = Lu(t) + Gf(t), 
u(0) = L~ l y 

with u(t) = L~ 1 y(t). Note that due to this transformation, we only obtain a convergence result 
for the original solution y(t) in the extrapolation space T>(L)*, i.e., with respect to the ||L _1 • || 
norm. 

First-order convergence of the exponential Lie splitting in D(L)* requires the regularity as- 
sumptions 

y^V^- 1 ), y'(0)=Ly + LGf(0)£V(L~<- 1 ), f G C 1+e ([0, T}; L 2 (d£l)) 
for some 7 > 0. Thus, the initial value yo and boundary data / need to satisfy 

(51) y eH 2 ^n), feC 1+9 ([0,T];L 2 (dn)) 
for some 7 > 0. 

Convergence of order 1.25 — e of the exponential Strang splitting in D(L)* requires 
y"(0) = L 2 y + L 2 Gf(0) + LG/'(0) G V(V~ l ), f G C 2+e ([0, T]; L 2 (dtl)) . 
Thus, the initial value yo needs to satisfy 

(52) y G V(L m ) n H 2+2 ^(n) = {<p G H 2+2 ^(n) ; ^\ dn = /(0)} 
for some 7 > 0. 

Considering the convergence in L 2 (il), i.e., measuring the error of the splitting methods ap- 
plied to (49) in the standard L 2 (Q) norm, both splitting methods suffer from an order reduction 
to j — e for all e > 0, as illustrated in Figure [4] below. This can easily be seen in the local error 
representation of the splitting methods, replacing the inhomogeneity g(t n ) with Lg(t n ) in (22) 
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and (38), respectively. For the Lie splitting and g replaced by Lg in the local error (22), the 
teTm.hr)C2g(t n ) can only be bounded by the loss of powers in h, as K^g{t n ) is only bounded for 
7 < j ~ e, i-e., 



h 2 \\)C2g(t n )\\ <h 1+ \- £ \\K,\ £ g(t r , 



Thus, the global error of the exponential Lie splitting in L?{£1) reduces to \ — e. Similarly, one 
can show the same for the exponential Strang splitting method. A numerical experiment where 
the inhomogeneity is simply chosen g(t) = 1 confirms this theoretical observation, see Figure |4j 
Again we summarize our findings in a corollary. 

Corollary 11. — Measuring the error with respect to the • || norm the exponential Lie 

and Strang splitting applied to (49) are convergent with order 1, respectively, 1.25 — e (for every 
e > 0) under the natural assumptions (51) and (52), respectively. With respect to the standard 
L 2 (f2) norm both methods suffer from an order reduction to 0.25 — e (for every e > 0). 



5. Numerical experiments 



In this section we numerically analyze the order of convergence of the exponential Lie (15), 



Strang (I24J) and Strang B (25 ) splitting to underline the theoretical convergence results derived in 
Section [3j The main interest lies in the sharpness of the convergence bounds for the exponential 
Strang splitting for inhomogeneous equations given in Theorem [7| as well as the convergence 
bound for the exponential Lie and Strang splitting for general inhomogeneous Dirichlet boundary 
conditions given in Corollary 11 In the numerical experiments we use standard symmetric finite 
differences for the spatial discretization of the split differential operators A(d x ) = d x ( y a(x,y)d x ^ 
and B(d y ) = d y (b(x,y)d b ). 



In Examples 12 and 13 presented below, we solve the inhomogeneous model problem (!]) with 
different choices of inhomogeneities ip. The numerical results confirm the error bound (26). In 



Example 14 we solve the inhomogeneous boundary problem (|49|) with the simplest choice of 

1. 



inhomogeneity g(t) 
Corollary [TT] 



The numerical result confirms the severe order reduction stated in 



Example 12 (Order reduction of the exponential Strang splitting) 

In this example we choose the following data in ([!]): 

C(d x ,d y ) = d x ((2xy + 3)d x ) +d y {(2xy A + l)d y ) 1 

8- 



wo{x,y) 



i i 

x(l— x) y(l—y) 



ijj(t, x, y) = x(l - x)y(l -y)+ te x y . 



Note that the regularity assumptions (12) and (14) are satisfied as in particular ^(0, -)\q^ = 0. 
Furthermore note that ip(t, •, -)\dn 7^ for t > 0. Thus, 7 + v must be chosen smaller than one 



quarter in (26) 



Hence, we expect first-order convergence of the exponential Lie splitting ( |15| ) and convergence 
of order 1.25 — e for the exponential Strang splitting method (24) applied to ([T]) formulated in 
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Figure 1. Numerical order of the exponential Lie, Strang and Strang B splitting applied 
to equation with the data as in Example 
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The error is measured in a discrete L 2 



norm. The slopes of the dashed and dash-dotted lines are one and 1.25, respectively. 
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Figure 2. Numerical order of the exponential Lie, Strang and Strang B splitting applied 
to equation with the data as in Example 12 The error is measured in a discrete V(L)* 
norm, i.e., with respect to \\L~ X • ||j,2. The slopes of the dashed and dash-dotted lines 
are one and two, respectively. 



L 2 (ft). On the contrary, we expect full-order convergence of both methods when formulating 
the evolution equation ([I]) in the extrapolation space (T>(L)*, \\L-~ 1 ■ ||), i.e., when the error is 
measured with respect to the ||L _1 • || norm, cf. Remark |)1 The numerical experiments go in line 
with our theoretical results, see Figure 111 for the L 2 (Q) and Figure^] for the T>(L)* formulation. 
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Figure 3. Numerical order of the exponential Lie, Strang and Strang B splitting applied 
to equation with the data as in Example 13 The error is measured in a discrete L? 
norm. The slopes of the dashed and dash-dotted lines are one and two, respectively. 



Example 13 (Full-order convergence of the exponential Strang splitting) 

Here we choose the following data in Q: 

C(d x ,dy) = d x ((2 X y + 3)d x )+d y ((2xy i + l)d y ), 



ip(t,x,y) 



i 



x(l—x) y(l—y) 



x(l - x)y(l - y)e t . 



Note that the regularity assumptions ( [12] ) and (14) are satisfied as tp(0, -, -)\dn = . In particular 
g(t) = ip(t, •) € 2?(L 1+7 ) for < 7 < - and all t. Hence, v can be chosen equal to one in (26). 



Thus, we expect first-order convergence of the exponential Lie splitting (15) and second-order 
convergence of the exponential Strang splitting method ( |24[ ) applied to (|lj). The numerical 
experiment goes in line with our theoretical results, see Figure [3j 

Example 14 (Inhomogeneous Dirichlet boundary conditions) 

Here we model the inhomogeneous boundary problem (49) with the data 



9(t) = l, y {x,y) 



i 



x(l—x) y(l—y) 



In Figure [4] we clearly see the severe order reduction of the exponential Lie and Strang splitting 
in L 2 (Q), which underlines the sharpness of Corollary 11 
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